Stripes and spin-incommensurabilities are favored by lattice anisotropics 
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Structural distortions in cuprate materials give a natural origin for anisotropies in electron prop- 
erties. We study a modified one-band t—J model in which we allow for different hoppings and 
antiferromagnetic couplings in the two spatial directions {tx ty and Jx Jy)- Incommensurate 
peaks in the spin structure factor show up only in the presence of a lattice anisotropy, whereas 
charge correlations, indicating enhanced fluctuations at incommensurate wave vectors, are almost 
unaffected with respect to the isotropic case. 
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The experimental finding of electron inhomogeneities, 
the so-called stripes, in doped Mott insulators has re- 
cently stimulated a great debate. The first experimen- 
tal detection of stripes in a family of La2-xSrxCu04 
(LSCO) cuprates was achieved in a Neodymium co-doped 
compound, Lai.6-xNdo.4SrxCu04. By using neutron 
diffraction measurements, Tranquada and co-workers, 
found that the magnetic peak at Q = (tt, tt) splits by 
a quantity x, giving rise to four incommensurate peaks 
around Q. Moreover, the Bragg peaks split by 2x around 
the F point. This fact has been interpreted with the 
formation of charge domain walls separated by antifer- 
romagnetic regions where no holes are present. More- 
over, the incommensurate magnetic peak reveals a tt- 
phase shift of the staggered magnetization when crossing 
the domain walls. Further studies Q show that similar 
low-energy magnetic peaks occur also in YBa2Cu306.6j 
indicating that the incommensurate spin peaks are com- 
mon features of all cuprates. The Nd-doping on LSCO 
produces a structural transition of the Cu02 plane from 
the low-temperature orthorhombic (LTO) to the low- 
temperature tetragonal (LTT) phase. ^] In the LTO 
phase the Oxygen atoms are displaced out the Copper 
plane and there is only one Cu— O bond length, see 
Fig. 0(a). In the LTT phase, instead, the Cu— O bonds 
are not equivalent in the two directions. In one direc- 
tion the Oxygen atoms are precisely in the Copper plane, 
while in the other the Oxygens are displaced out of the 
Copper plane, see Fig. |l|(b). The Cu— Cu hopping de- 
pends on the Cu— O bond and it is isotropic in the LTO 
phase and anisotropic in the LTT one. 

In this letter, we show that the anisotropies in the 
Cu02 may help in determining incommensurate electron 
correlations, j^,^ In particular, the anisotropy in the 
antiferromagnetic coupling favors huge incommensurate 
spin peaks. These peaks, although reduced, persist also 
when considering the anisotropy in the hopping, indicat- 
ing that the structural transition may be responsible for 
the electron incommensurability, even if other physical 
effects may cooperate in the occurrence of the striped 



phase. The fact that no sign of enhanced charge fluctua- 
tions are detected in the presence of spatial anisotropies 
does not contradict the experimental finding, which are 
sensible to the spin degrees of freedom and not to the elec- 
tronic charge. Finally, we show that the superconducting 
pairing is suppressed only along the stripe, supporting a 
coexistence of incommensurate spin fluctuations and su- 
perconductivity, as found in Ref. [^. 

We consider a two-dimensional t—J model with differ- 
ent coupling in the x and y directions: 

i,fi—x,y 

- Yl tf,dl^c^+^,^„ + H.c, (1) 
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where ^ = c| g. (1 — n^^^), and Si are the electron 
density and spin on site i, respectively. The anisotropies 
are given by the fact that — a^t and = I3^J . 

We use different Quantum Monte Carlo (QMC) tech- 
niques to estimate the ground-state properties of the 
Hamiltonian ([^). In particular, we apply the recent 
QMC method based on the application of p Lanczos steps 
to a given variational wave function. Q Moreover, we 
also used the fixed-node j|] and stochastic reconfigura- 
tion (SR) approximations. We consider the p = 
and the p = 1 state as the guiding wave function for 
the fixed-node method, and, hereafter, the symbols FN 
and FNLS indicate the fixed-node approximation with 
the p — and p = 1 wave functions, respectively. Al- 
though these techniques provide a considerable improve- 
ment of the simple variational calculation, they crucially 
depend on the choice of the guiding wave function. Both 
the fixed-node and the SR are exact both in the limit of 
strong anisotropy {t^ — Jx = and any doping) and in 
the low-doping limit (any anisotropy) and therefore they 
represent reliable approximations for this problem. 

The best variational wave function is the projected 
BCS state 
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where Vn is the projector onto the subspace of N parti- 
cles, Vg is the GutzwiUer projector, which forbids dou- 
bly occupied sites, J = exp(^j ^ Wij/ii/ij) is a Jas- 
trow factor, defined in term of the hole density at site 
i hi = (1 — ni|)(l — rii^i), and Vij are variational param- 
eters. The variational parameters fij represent the pair 
amplitude of the BCS wave function. 
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FIG. 1. A schematic representation of the Cu02 plane 
in the LTO (a) and in the LTT (b) lattice structure. The 
black circles represent Coppers, which are in the plane in 
both structures. Large circles represent Oxygens, the symbol 
+ (— ) indicates that the atom is shifted above (below) the 
Copper plane. 

The wave function (^) naturally describes a Resonat- 
ing Valence Bond (RVB) state in which preformed pairs 
become superconducting by doping. |^,|l^ Whenever the 
lattice does not break any translational and rotational 
symmetries and J ^ 0.5t, the wave function is an 
exceptionally good approximation of the ground state of 
the t—J model, j?) For instance, in the case of 4 holes on 
a 26 lattice and J = 0.5t, the overlap of this wave func- 
tion with the exact ground state is 0.88. In the isotropic 
model, by using QMC and starting from this very accu- 
rate wave function, it is then possible to obtain almost ex- 
act ground-state properties also for large lattices, indicat- 
ing a small but finite superconducting order parameter 
and no charge inhomogeneities. [pl| On the other hand, 
density matrix renormalization group (DMRG) calcula- 
tions for J ^ 0.5t are interpreted as a strong tendency of 
the ground state toward the formation of static stripes. 

These results, however, were obtained for rectangular 
clusters and a particular choice of the boundary condi- 
tions: periodic (open) along the short (long) direction. 
The first choice breaks the rotational symmetry, the sec- 
ond one breaks also the translational invariance. 

In order to understand how stripes are stabilized, it is 
important to release in the model only the rotational in- 
variance, without spoiling the translational symmetry. In 
fact the LTO — > LTT transition leads to a conformation 
which naturally breaks the rotational invariance, preserv- 



ing the translational symmetry of the underlying lattice. 
We begin by considering a rectangular lattice 4x8 with 
periodic boundary conditions (PBC) on both directions 
and 4 holes, for the case with no anisotropy in the cou- 
plings, = = 1 and J = OAt. In Fig. ||(a,b), we 



report the density-density correlations N{q) 



for the variational wave function with p = 0, for the FN, 
FNLS and SR approximations. In this case the varia- 
tional state does not present any particular structure in 
the density correlations, and the p — 1 and 2 results 
do not change the overall nature of the state. Instead, 
the more involved FN, FNLS and SR approaches give 
rise to a huge peak at Qc = (0,7r/2). The FN, FNLS 
and SR methods are able to detect the right tendency of 
the charge arrangement, because the approximate ground 
state sampled by these techniques can allow the forma- 
tion of many-holes bound states, not present in the wave 
function (0), containing only two-body correlations. 




1 

^ 0.8 
^' 0.6 

0.2 








-_ (c) 









M 



M 



FIC. 2. N{q) and S{q) for 4 holes on a 4 x 8 lattice and 
J = OAt. (a) Variational with p = 0, (b) FN (empty squares), 
FNLS (full squares) and SR (full circles) results for N{q). (c) 
Variational with p — 0, (d) FN (empty squares), FNLS (full 
squares) and SR (full circles) results for S{q). F — (0,0), 
M = (^,^),V = (^,0), V = (0,^). 

From the experimental point of view, the most impor- 
tant signature of stripes is the appearance of well-defined 
incommensurate peaks in the magnetic structure factor 
S{q^uj) at small energies. It is really impressive that, 
by using the FN, FNLS and SR approaches, we find the 
same fingerprint of stripes in the equal-time correlations 
S{q) = (SgStg), as reported in Fig. ||(c,d). For the most 
accurate calculations, the Q-peak splits into two peaks 
at Qs = (tt, TT ± 7r/4), indicating that the modulation in 
the spins is twice the one in the charge. 

This calculation clearly shows that stripes can be de- 
rived by weak perturbations over the state of Eq. (||): 
many pairs may bound together along a preferred direc- 
tion, and are stabilized in the t—J model, if the rotational 
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symmetry of the Cu02 plane is explicitly broken. 
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FIG. 3. ^(si) for 8 holes on a 8 x 8 isotropic lattice and 
J = 0.4t. Variational with p = (empty triangles), FN (full 
triangles) and SR (full circles) results are reported. In the 
inset: the same for S{q). 

In the following we compare square lattices with PBC 
on both directions with and without anisotropics. In 
the latter case, for moderate hole doping ^ ~ 0.1 and 
J ~ 0.5t, the wave function (|^) is an accurate approxima- 
tion of the ground state: as shown in Fig. |^ there are no 
changes in the correlation functions by improving the ap- 
proximation used. As emphasized in Figs. |^ and^, within 
our approach, there is a strong influence of the bound- 
ary conditions for the stabilization of stripes. Within 
the DMRG, it was not possible to realize the important 
role played by the boundary conditions. Indeed, it is 
well known that DMRG calculations are not accurate for 
PBC, especially on two-dimensional clusters, where it is 
not possible to reproduce the homogeneous ground state 
even on the 6x6. Although for few chains our approach 
is less accurate - but qualitatively correct - than DMRG, 
Ip^ it is more suitable for two-dimensional lattices. The 
fact that QMC is competitive with DMRG is confirmed 
by the fact that on an 8 x 8 lattice with open boundary 
conditions (OBC) on both directions (where the DMRG 
is much more accurate than PBC), J = 0.4t and 8 holes, 
the SR variational energy, E = — 39.296(6)^, is lower than 
a "state of the art" (4200 states kept) DMRG calculation, 
E — — 39.2503i, and very close to the extrapolated result, 
consistent for both techniques, E = — 39.44(4)^. |l5| 

The presence of anisotropics strongly affects the out- 
come given by Fig. ||. Indeed, by comparing Fig. ^ with 
the analogous one without anisotropy (Fig. ||) , it appears 
that the spin-spin correlations are strongly affected by 
the anisotropy: though at the pure p = variational 
level S{q) has a broad peak around the antiferromag- 
netic vector Q, within the SR technique, incommensu- 
rate peaks at (tt, tt ± 7r/4) show up, see Fig. ^(a,b). We 
expect that the exact value of these peaks is underesti- 
mated and that they are much more pronounced in the 



exact ground state. Indeed, whenever the ground-state 
correlations are qualitatively different from the ones of 
the projected BCS wave function, by using the SR (or 
the fixed-node) approach it is possible to detect the most 
relevant changes in the correlations. The incommensu- 
rate peaks for 8 holes on 8 x 8 and 18 holes on 12 x 12 
lattices are consistent with two half-filled and 3/4-filled 
stripes, respectively. It is worth noting that 3 half-filled 
stripes in the 12 x 12 give rise to 3 7r-shifts which are not 
compatible with PBC. Thus, it is not clear at present 
what the is most stable stripe filling in the thermody- 
namical limit (from half to completely filled stripes |^ ) . 
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FIG. 4. (a): S(q) for = 1.1, Oy = 0.9 and = 
(case 1). Empty and full circles: variational with p — and 
SR results for 8 holes on a 8 x 8 lattice. Full triangles: SR 
results for 18 holes on a 12 x 12 lattice, (b) : S(q) for 8 holes on 
the 8x8 lattice with — 1, f3x = 1.2 and f3y — 0.8 (case 2). 
(c): the same as (a) for N{q). (d): the same as (b) for N{q). 
In the inset: the energy per site extrapolation for p = 0, 1, 2 
Lanczos steps for the case 1 (continuous line) and for the case 
2 (dashed line), for the 8x8 lattice. The extrapolated points 
are also shown. 

The effect of the lattice anisotropy on the charge 
correlations is much less evident and N{q) does not 
show sizable differences between different approaches, see 
Fig. ^(c,d). We stress that the neutron scattering is only 
sensible to the electronic spin and there is no direct infor- 
mation on electronic charge. The finite peak at small-^/ in 
the N{q) may be attributed to dynamical fluctuations of 
the stripes, which do not affect the S{q) incommensurate 
peak, but destroy a coherent response in the static charge 
correlations. In order to have this outcome, it is impor- 
tant that the spins across the stripes are strongly antifer- 
romagnetically correlated. ||l^ We have verified that, in 
the anisotropic case, the spin-spin correlation across an 
hole has a strong antiferromagnetic character. The above 
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scenario is confirmed by the fact tliat the anisotropy in t 
favors the transverse motion of the domain walls, yield- 
ing a suppression of the incommensurate peak of S{q), 
as shown in Fig. §(a,b) for the 8x8 case. Well defined 
static stripes are instead defined in the system by con- 
sidering OBC. As shown in Fig. ||, in the 8x8 cluster 
there is a clear tendency to formation half-filled stripes, 
although in the variational calculation no signature of 
inhomogeneities is present. For 4 holes there is a single 
stripe at the center of the lattice, whereas for 8 holes, two 
of such features clearly appear, indicating the half-filled 
character. 
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FIG. 6. (a): pairing correlation function Aij^k,i obtained 
with SR for the isotropic 8x8 lattice of Fig. ^ (circles) and the 
anisotropic one of Fig. ^b,d) (triangles and squares), (b): the 
same for the 4x8 lattice of Fig. ^. The pairs {k, I) (thin lines) 
are obtained by moving the pair {i,j) (thick line) parallel to 
the X or the y axis, d is the distance between pairs. The 
direction of the stripes is defined by the arrow. 
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Rice, M. Calandra, A. Parola, M. Fabrizio, A. Bianconi 
and C. Morals Smith for useful comments. This work has 
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FIG. 5. On the left: average hole density profile for 4 holes 
in a 8 X 8 lattice with OBC and J = OAt, = 1, = 1.2 
and Py = 0.8 (case 2 of Fig. ^), for the p — calculation 
(upper- left panel) and for the SR technique (lower-left panel) . 
On the right: the same for 8 holes. 

In order to study the role of the anisotropics on the 
superconductivity, we consider the pair-pair correlation 
function Aij^k,i, which creates a singlet in the sites i 
and j and destroys it in the sites k and I. The results 
are reported in Fig. ^. The effect of the anisotropy is to 
depress the pairing function along the stripe, whereas the 
pairing remains almost unchanged in the perpendicular 
direction. This effect is particularly strong in the 4x8 
lattice where pairing correlations are suppressed by more 
than an order of magnitude along the stripe. 

In conclusions, in the presence of anisotropies, our find- 
ing is consistent with fluctuating stripes, where the tt- 
phase shift gives rise to incommensurate peaks in the 
S{q) and to an almost featureless N{q). We finally re- 
mark that stripes do not necessarily suppress supercon- 
ductivity, but large pairing correlations can be obtained 
in the direction perpendicular to the stripes. On the 
other hand, d-wave superconductivity is a very robust 
property of an isotropic doped antiferromagnet, implying 
that the RVB scenario represents a reasonable explana- 
tion of high-temperature superconductivity. 
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